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Abstract 

We present a new topological method for the study of the dynamics 
of dissipative PDE's. The method is based on the concept of the self- 
consistent apriori bounds, which allows to justify rigorously the Galerkin 
projection. As a result we obtain a low-dimensional system of ODE's 
subject to rigorously controlled small perturbation from the neglected 
modes. To this ODE's we apply the Conley index to obtain information 
about the dynamics of the PDE under consideration. 

We applied the method to the Kuramoto-Sivashinsky equation 

lit = {u 2 )x — Mm — VUxxxx, u(x, t) — u{x + 2ty, t), u(x, t) = — u{— X, t) 

We obtained a computer assisted proof the existence of the large num- 
ber fixed points for various values of v > 0. 
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Table 1: Coefficients for the function u(x). 



1 Introduction 

Even in the setting of infinite dimensional dynamics many of the dynamical 
objects of interest arc low dimensional, e.g. equilibria, periodic orbits, connect- 
ing orbits, horseshoes, etc. In this paper we introduce techniques which, in 
principle, allow for the rigorous verification of such solutions for a wide variety 
of partial differential equations. Our approach is to combine rigorous computer 
calculations with topological invariants to obtain accurate existence statements. 
To demonstrate these techniques we have chosen to study the the Kuramoto- 
Sivashinsky equation jjjl |l9| 

tit = -vu xxxx - u xx + 2uu x (t, x) £ [0, oo) X (— 7T, 7r) (1) 

subject to periodic and odd boundary conditions 

u(t, — tt) = u(t, 7r) and u(t, — x) = — u(t, x). (2) 

The following theorem is a prototype for the results which can be obtained. 

Theorem 1.1 Let u{x) = X)fc=i a fe sin(fcx) where the are given in TableQ 
Then, for v = 0.1 there exists an equilibrium u*(x) for (^) such that 

\\u* - ii|| L 2 < 2.71547 x 10~ 13 and | \u* - u\ \ c o < 2.06706 x 10~ 13 . 

Having stated this theorem we now try to put the result into the context 
of the goals of our methods. To begin with it needs to be emphasized that the 
computations which lead to this result are rigorous in the sense that we have 
employed interval arithmetic to overcome all errors introduced by the fact that 
we are using floating point arithmetic in our calculations. 

As it will become clear in the later sections, this result is obtained by study- 
ing the full partial differential equation rather than attempting to solve a bound- 
ary value problem. While from the point of view of traditional numerical anal- 
ysis this approach may appear inefficient, it is an important point. To be more 
precise, our method does not attempt to directly approximate any particular 
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Stationary Solution to Kuramoto-Shivashinsky v = 0.1 



Domain is [- 71,71] 



Figure 1: The equilibrium solution u*(x) for v = 0.1 



solution to the partial differential equation. Rather we essentially compute the 
Conley index of a compact region, called an isolating neighborhood, of phase 
space. The diameter of this region provides the error bounds stated in the 
theorem. The index guarantees the existence of the equilibrium solution. 

The Conley index theory is a far reaching topological generalization of Morse 
theory. In particular, this index can be used to prove the existence of periodic 
orbits, connecting orbits, and chaotic dynamics ||, |2^, O, |^, It has been 
numerically observed that for various parameter values (jjl) contains these types 
of dynamical objects. In principle, combining earlier rigorous numerical methods 
| jl4[ 15, [l?], [l], 0, [24|, [2^] with the techniques described in this paper and the 
above mentioned index theorems will lead to rigorous proofs of the existence of 
periodic orbits and even chaotic dynamics. However, we do not pursue these 
more complicated structures in this paper for two reasons. First, finding the 
appropriate isolating neighborhoods is more complicated in these cases and our 
goal here is to emphasize the fundamental ideas associated with the methods. 
The second, and more important point, is that a straightforward application 
of the earlier numerical methods would lead to large computations - which we 
believe can be avoided by alternative methods (see for example p3| ). This latter 
point is currently being investigated. 

Returning to our discussion of Theorem |l.l| , an obvious question concerns the 
stability of u* . For this we have no definitive answer. As was indicated before 
our method docs not directly approximate u* and therefore we do not obtain 
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uniqueness results or hyperbolicity results. On the other hand, being a gener- 
alization of the Morse index the Conley index does contain some information 
about the stability or instability of the dynamics in the isolating neighborhood. 
Thus, what can be asserted is the following. Assume that u* is a hyperbolic 
fixed point, i.e. all eigenvalues have nonzero real part, and that u* is the only 
solution which remains within either the L 2 or C° bounds of u for all time, 
then u* has exactly two unstable eigenvalues, i.e. its unstable manifold is two 
dimensional. We hope to treat this problem in a subsequent paper. 

It should be mentioned that even though we are doing the computations 
via an approximation of the full partial differential equation, we never integrate 
the equations. Rather, as will be made clear in Section 2 the computations are 
reduced to solving a set of inequalities. It is for this reason that we are able to 
get such sharp bounds on the equilibria. As the following theorem demonstrates 
we can, in fact obtain bounds on the level of the floating point accuracy. 

Theorem 1.2 One can compute a sequence a\, a^, . . . , 030 and the function 
u(x) — X)fe=i °fc sin(fcx) ; such that for v = 0.75 there exists an equilibrium u*(x) 
for (@j such that\\u*~u\\ L 2 < 1.26281 x 10" 15 and ||u*-u|| c o < 9.57396x 1CT 16 . 

By now it is a well demonstrated principle that the asymptotic behavior of 
a wide variety of infinite dimensional dynamical systems is finite dimensional 
|j| ||, |2lJ . The Kuramoto-Sivashinsky equation (|l|) is a particularly well studied 
example of such a system 0. In fact, it is known that (|l|) posses an inertial man- 
ifold and therefore, that there exists a family of ordinary differential equations 
that exactly describes the asymptotic dynamics. Unfortunately, the estimates 
for the dimension of these manifolds make them impractical for our purposes 

We mention these methods to emphasize that our approach does not di- 
rectly make use of any of these results. What appears to be essential for our 
techniques is that the spectrum of the linear operator for the evolution equation 
is not clustered near the imaginary axis. This is in contrast to the inertial man- 
ifold techniques which strongly rely on gap conditions of the spectrum or cone 
conditions from the flow. Our approach is to use the computer to restrict our 
attention to that portion of phase space in which the desired dynamics (for this 
paper the fixed points) occur. Obviously, by restricting the phase space one can 
get much better estimates. This sets up a loop by which one can continuously 
improve the estimates until the desired bounds are reached. 

Our analysis of the fixed points for (|l|) was motivated in part by the work 
of Jolly, Keverkidis, and Titi pj. In particular, using a 12 mode traditional 
Galerkin approximation, they produced a bifurcation diagram for v £ (0.057, 00). 
We used their reported solutions to test our methods. In particular, as is indi- 
cated below we were able to find and prove the existence of an equilibrium point 
on each of their stable branches. Unfortunately, we used a fairly primitive search 
procedure and therefore missed a few unstable branches. Our expectation is that 
by combining our methods with a continuation package, one could produce a 
rigorous bifurcation diagram with fairly precise bounds in a computationally 
inexpensive manner. 



4 



Below we include some of the steady states we found. 



• v = 0.5. Two stable unimodal fixed points 

• v = 0.3. Two stable unimodal fixed points 

• v = 0.i27, v = 0.1 25. A stable and unstable bimodal fixed point. Negative 
branch is stable, positive one is unstable with apparently two-dimensional 
unstable manifold. 

Our primitive search procedure did not find a solution on the bi-tri branch. 

• v = 0.1 An bimodal stable and unstable (2 unstable directions) and two 
unstable trimodal fixed points (both with 1-dimensional unstable mani- 
fold) 

We did not find an unstable branch connecting bi-tri branch with quadri- 
modal branch. 

• v — 0.08 A bimodal stable (neg. branch) and unstable (2 unstable direc- 
tions) fixed points. A pair of stable fixed points close to R^ti (see p0[). 
A pair of unstable trimodal fixed points (1 unstable direction). 

We did not find a branch connecting bi-tri and quadimodal branches. 

• v = 0.0666.., v = 0.063 Two unstable bimodal points, two stable trimodal 
points and two stable solutions apparently belonging to the giant branch. 



We are lacking two unstable branches which are present in 1 10 



• v = 0.062, Two stable trimodal points and two stable points from giant 
branch 

• v = 0.045, Two stable points from giant branch and pairs of unstable tri- 
and quadrimodal fixed points 

• v = 0.04, Two stable giant fixed points. Two stable quadrimodal fixed 
points. Two unstable trimodal points. 

• v = 0.029 Two unstable quadrimodal points. 



2 The Method 

Our method begins with the reduction of the full dynamical system to a lower 
dimensional system which can be studied numerical. In particular, we begin 
with a nonlinear evolution equation in a Hilbert space H (L 2 in our treatment 
of Kuramoto-Sivashinsky) of the form 

where domain of F is dense in H. Furthermore, we assume that | i = 
0, 1, . . .} forms a complete orthogonal basis for H . 
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In the case of the Kuramoto-Sivashinsky equation F(u) = Lu + B(u,u), 
where L is a linear part and B is a nonlinear part, the functions {ifi} are chosen 
to be eigenvalues of L. 

Fix to £ N and let 

P = P m : H — » A m = X 

be the orthogonal projection from H onto the finite dimensional subspace span- 
ned by {0i, 02, ■ ■ ■ , 0m}- Let 

Q = Q m ■■= i - p ■ h y = y m 

be the complementary orthogonal projection. Finally, let 

A k : H R 

be the orthogonal projection onto the subspace generated by 0&. 

Given u £ H, let Pit = p and Qu = q. Equation (S) can be rewritten as 

f = PF(p,q) (4) 
§ = Q^(P,9) (5) 

The strategy adopted is fairly simple: study the dynamics of the low di- 
mensional Galerkin projection (^) to draw conclusions about the dynamics of 
(§). Before tur ning to the precise conditions, consider the following heuristic 
description of our approach. 

Let W C X = X m . For j > to, let Wj C X 3 such that P(p- 1 (W j )) = W, 
(i.e. Wj = W © (I - P)Wj). Similarly, let V C Y and set Vj = Qj(V). 
Furthermore, given qj £ Vj assume that limj^oo \ \qj\ \ — 0. Our only knowledge 
concerning the higher order modes or "tails" of the solutions is that they project 
into V. This implies that our knowledge of the vector field is reduced to the 
following differential inclusion 

^£PF(p,V) 

where p £ W. Numerical calculations on this equation are used to find topo- 
logical invariants (the Conley index, the fixed point index) which guarantee the 
existence of specific dynamics, e.g. fixed points, periodic orbits, symbolic dy- 
namics, positive entropy, etc. It is simultaneously argued that the topological 
invariant is the same for any Galerkin system of the form 



where pj £ Wj. Thus, the same dynamical object exists for each sufficiently 
high Galerkin approximation. Finally, it is shown that the limit of these objects 
leads to the desired dynamics for the full system (|J). 
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2.1 Self-consistent Bounds 



As one might expect the orthomormal basis {4>i} and the sets W and V must be 
chosen with care. The first issue that needs to be dealt with is analytic in nature 
- solutions to the ordinary differential equations must approximate solutions of 
the partial differential equation. This leads to the following definition. 

Definition 2.1 Let m,M G N with m < M. A compact set W C X m and 
a sequence of pairs {a^ G K | < a£, k G N} form self-consistent apriori 
bounds for (|^) if the following conditions are satisfied: 

CI For k> M, a~ < < a+. 

C2 Let afc := max|a^| and set u — X^feLo a k4>k- Then, u G H. In particular, 
\\u\\ < oo. 

C3 The function u h- ► -F(it) is continuous on 

oo 

fc= 771+1 

In practice C rifeLi[ a fc ' a fc]* Given self-consistent apriori bounds W and 
K !• let 

OO 

v~ n K,a+]cy m . 

/c— m+l 

Our goal is to numerically solve (^) on and draw conclusions about the 
dynamics of (||) on the set W ® V C H. To do this we will make use of the 
following results, the first two of them are obvious. 

Lemma 2.2 Given self- consistent apriori bounds W and {a^ }, W © V is a 
compact subset of H. 

Lemma 2.3 Given self-consistent apriori bounds W and {a^ }, W © V, then 
lim Q n (F(u)) = 0, uniformly for u G W © V 



Proposition 2.4 Let W and {a^ } be self-consistent bounds for M). A function 
a : [0, T] -> W © V groen &y 

OO 

fe=0 

is a solution to if and only if for each k G N and < G [0, T] 

^ = A fe F(a). (6) 
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Proof. (=>) This direction follows directly from the projection of (g) onto 
each of the basis elements. 

(<=) Assume that (pf) is satisfied for each k S N and all t 6 [0, T]. Let 



First observe that from C3 it follows immediately that J2T=i MT'fik ~ ^( a ) ^ 
H. 

It needs to be shown that 



da a(t + h) — a(t) 

— = lim 

at h^o h ^-^ at 

fc=0 



Eda k 



This is equivalent to showing that 



lim 

h-*0 



1 00 00 d 

- ^(a fe (< + h)- a k {t))4>k 



k=l 



k=l 



= 



for all i e [0,T]. 

Fix /i > 0, then for any n E N 



! x; oo 

- >J(ofc(t + h)- a k (t))4>h - 



fe=X 



fc=i 



< 



- n n 

- 5^(o k (* + h)- a k (t))cf> k - J2 
fe=l fc=i 



fc=n+l 



dak 
~~dt 

Eda k 

k=n+l 



<Pk 



We will estimate the three terms on the right hand side separately. From 
lemma 2.3 it follows for a given e > there exists no such that n > no implies 



E 

k=n+l 



dak 
~~dT 



= \Q n (F(a))\ < e/3. 



From now on fix n > n . Again lemma 2.3 and the mean value theorem implies 
^ j{a k {t + h) - a k {t))(j) k 



k=n+l 



E + 



eft 

k=n+l 

\Q n (F(a(t + 6 k h))\ < e/3. 



Finally, for h sufficiently small, 

1 n n d 

- £> fc (t + fc) - a k (t))J> k -Y,-JT 



k=l 



k=l 



< e/3 



and hence the desired limit is obtained. 
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2.2 Conley Index 

Proposition |2.4| indicates that given self-consistent apriori bounds W and {a^}, 
finite time solutions to @ are solutions to the full partial differential equation. 
Thus, the goal of this paper is to find solutions to (||). Of course, numerically one 
can only study (||) restricted to W and then argue that the resulting numerical 
solution is an approximation to a solution to (^). Hence, rather than attempting 
to approximate specific trajectories in W directly, the objective is to compute a 
Conley index for (Q) and then show that this index information is sufficient to 
guarantee a solution for (||). 

In order to describe this index the following definitions are needed. Let 
p:Ix IR™ — > IR™ be a continuous flow generated by a differential equation 
z = f(z). 

Definition 2.5 A compact set N C R™ is an isolating neighborhood if 

Inv (N, ip) :={zeN \ <p(R, z) C N} C hit N. 

If, in addition, for any z £ ON, there exists t z > such that 

tp((0,t z ),z)nN = $ or <p((-t z ,0),z)nN = <b, (7) 

then N is an isolating block. Given an isolating neighborhood N, the associated 
maximal invariant set Inv (TV, ip) is an isolated invariant set. 

The easiest way to verify the existence of an isolating block is through local 
sections. 

Definition 2.6 EE C is a local section for ip if for some e > 

<p:(-e,e)xE^tp{(-e,e),S) (8) 

is a homeomorphism and <p({— e, e), 2) is an open subset of M" . 

A special form of local section is a hypersurface which is transverse to the 
flow. More formally, let S C M be an n — I dimensional manifold with normal 
vector ii(z) at z G S. Sisa local section if for each zeS, 

H(z)-f(z)?0. (9) 

It is straightforward to check that N is an isolating block if dN can be 
written as the union of the closure of local sections with the property that (0) 
is satisfied at every point in the intersection of the closure of the sections. 

In this paper the focus is both on proving the existence of equilibria and 
providing tight bounds on the location of the equilibria. To do this requires 
have good isolating blocks. With this in mind consider the linear ordinary 
differential equation 

z = Bz, ze R n . (10) 
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Assume that the origin is a hyperbolic fixed point. Without loss of generality 
it can be assumed that B is in Jordan normal form. Generically, to each real 
eigenvalue there is associated a 1-dimensional eigenspace and to each pair of 
complex conjugate eigenvalues there is an associated 2-dimensional eigenspace. 
Thus, IR™ can be decomposed into the product of eigenspaces, i.e. 

R n = Vi X V 2 X • • • X V k 

where Vi is either IR or IR 2 . In what follows we will use the following notation, 
Zi G Vi, i = 1, . . . , k, and if Vi = R , then Zi = (xi,yi). 

Our interest is not on the dynamics of (0) on the entire phase space, but 
rather on a prescribed compact subset. Since our goal is to understand the 
equilibria of ( |To| ) consider a neighborhood of the origin, 

N = h x h x ■■■ x h 

where 

j . = f KM K < o <bj if Vi = R, 

1 ' I {(x it yi) S V t | y f xf+v! < h, h > 0} if Vi S R 2 . 

The following result is obvious, but to make a point crucial to the results of 
this paper we will provide the proof. 

Lemma 2.7 The compact set N is an isolating block for 

Proof. Since B is in Jordan normal form the system decouples according to 
the decomposition R" = V\ x V 2 X ■ ■ ■ x Vk- 

If Vi = M, then ( |Toj ) reduces to Zi = \zi. Since B is hyperbolic ^ 0, and 
hence at Zi — bf (0) becomes 

Kbf + 0. 

If Vi = R 2 , then © reduces to 

Vt = -PiXi + aiyi 
where by hyperbolicity cti ^ 0. So again for \fx\ + yf = bi (^|) becomes 
(xi, yi) ■ (ctiXi + 0iyi, -PiXi + aji/i)* = ^ 0. 

I 

To see why this trivial argument is of importance, consider the more inter- 
esting example of 

z = Bz + f(z) + E(z) (11) 

where / : IR" — ► IR™ is o(||z|| 2 ) at and E represents a known bounded error. 
In our situation E arises from numerical errors and approximations. More 
precisely, we assume that there are known constants Ci such that 

SUp ||JSj(2)|| < Ci. 
z£N 
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Observe that a sufficient condition for TV to be an isolating block for (|ll|) is 
the following: for each i such that Vi = M, 

Xibf + fi(z) + Ei(z) (12) 

has the same sign as Xibf over the set {z G TV | Z{ = bf}, and for each i such 
that Vi = M 2 , 

(xi,yi) ■ (otXi + Pyi + fn(x) + E il (x),~Px i + ay { + fi 2 (x) + E i2 (x)Y (13) 

has the same sign as Oj over the set {z G TV \ yj ' x\ + yf — b t }. 

For the linear case the eigenvalues of B are assumed to be known and the 
bf can be chosen arbitrarily. Therefore, one can also interpret ( |l2"| ) and ( p"3| ) 
as providing a set of inequalities that if simultaneously solved for bf provide 
an isolating block even in the context of numerical errors and approximations. 
In particular, finding isolating blocks need not involve numerically solving the 
ordinary differential equation. 

In itself the knowledge that TV is an isolating block does not imply anything 
about Inv (TV, ip) . To gain information concerning the isolated invariant set we 
will make use of the Conley index. For our purposes we need only a very small 
portion of the index theory and so we give a minimal operational definition (see 
[g f2Cj |l8[ (| H for further information). 

Definition 2.8 Let TV be an isolating block and let dN = L + U L~ where L 
are closed sets. Furthermore, assume that z € L~ implies that 

y((O,e),z)nTV = 

for a sufficiently small e > 0. Similarly, assume that if z G nL + , then 

<^((-6,O),z)nTV = 
for a sufficiently small e > 0. The Conley index of S = Inv (TV, ip) is 

CH*(S) := H*(N,L). 

No knowledge of relative homology groups is required for the applications 
described in this paper. The following theorem gives a formula for the index of 
a hyperbolic fixed point. 

Proposition 2.9 Let q be the number of eigenvalues of B with positive real 
part. Assume that for all i = 1, . . . , k either the condition associated with |IJ 
or the condition associated with fildj ) are satisfied . Then, 

Z if j = q 
otherwise. 



CHj{lTiv(N,ip)) 



The following theorem due to McCor d [p"3"[ Corollary 5.9], indicates that if 
the Conley index is a that of Proposition 2.9, then there exists a fixed point in 
TV. 
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Theorem 2.10 If the Conley index has the form 



CHj(Inv(N,<p)) = 




for some q, then N contains a fixed point. 

To indicate how these index ideas will be used in this paper let us return to 
the system ([H]). Observe that the only assumption on the error term E was 
that it is bounded, therefore, it is no longer apriori true that the origin is a fixed 
point or that there even exists a fixed point to ([ll]) . On the other hand the sets 
N and L remain unchanged. Therefore, the Conley index implies the existence 
of the fixed point. 

2.3 A Singular Perturbation Result 

As was indicated in the previous section, it is possible to find an isolating block 
for a finite dimensional ordinary differential equation about a fixed point by 
solving an appropriate set of inequalities. However, to do this requires a good 
estimate of the location of the fixed point, knowledge of the eigenvalues, the abil- 
ity to evaluate the nonlinear terms, and estimates of associated errors. There- 
fore, the dimension to which one can hope to apply this procedure is obviously 
limited. In this section we will describe a singular perturbation result which 
allows one to "lift" the index computations of the previous sections to arbitrary 
dimensions. 

The definition of self-consistent bounds related individual solutions of the 
infinite family of ordinary differential equations to solutions in partial differential 
equation. We now need to extend this definition in order to know that the 
index computations we perform for the finite dimensional approximation have 
implications for the partial differential equation. 

Definition 2.11 Let m, M £ N with m < M. A pair of compact sets N C 
W C X m and a sequence of pairs {a k € M | < at, k € Ft} are topologically 
self-consistent if W and {a^} are self-consistent apriori bounds and the following 
conditions are satisfied. 



For Kuramoto-Sivashinsky we will make use of the following stricter form of 



C4 Let u S W © Uk=m+i[ a k ' a tl Then ' for k>m 

A k u = <zj =*> A k F{u) j= 0. 

C5 N is an isolating block for (Q) for all q e Ilfe>m[ a fc > a t 



(14) 



C4. 



C4a Let u e W © ft 



OO 



[a k , at]. Then, for k > m 



k=m-\-l 



A k u = a+ A k F{u)) < 

A k u = a^ A k F(u)) > 0. 



(15) 
(16) 
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Using the line of reasoning that was described in the analysis of (|TT|) , con- 
dition C5 can be replaced by the following assumption. 

C5a Let N be an isolating block for (Q). Let v ± {p) be the outward normal at 
p £ L ± . If u G W © ]lfclm+i[ a fe i a fc 1 such that Pu£L ± , then 

PF(tt) ■ f + (Ptt) < PF(u)-u-(Pu) >0. 

We shall now discuss two singular perturbation results. The first allows one 
to lift isolating blocks. 

Theorem 2.12 Let m, M £ N wifA m < M. A ssume N £ W £ X m and 
the sequence of pairs {a^ £ R | a t < a^ , fe € N} are topologically self- 
consistent. Fix an integer r > m. Then for any q = J2T=r+i ( lk<t l k, such that 
q £ n^ r+1 [a^, ajj"] and q k = for k > M the set 

N :~ N x [a m+1 , a+ +1 ] x [a m+2 , a+ +2 ] x ■ • • x [a", a+] 

is an isolating block for the system of equations 

k 

x k = AkF(}^ Xj(j)i + q) fc=l,...,r (17) 

i=l 

where x £ R r . 

Proof. Let u = (w,v) £ W © Ilfc=m+i[ a fc i 4\- From CI it follows that 
u + q £ W © TlT= m +ii a k > a t]- If w G 9iV, tnen either w is in cW or v is in 
^rifc=m+i[ a fe > a t ]• ^ n tne nrst case ^5 forces the vector field to be transverse 
at the boundary. In the second case transversality follows from C4. I 

Remark 2.13 For r > M equations ftT\ ) are the Galerkin projection of u = 
P(u). 

The direction of the vector field influences the index computation. With this 
in mind define 



dir(fc) := 



'-1 


if A^u = 


4 - 


> A k F{u) 


< and 




A k u = 


a k -- 


» A k F{u) 


) > o 





if AkU — 


4 - 


> A k F{u) 


< and 




A k u = 


a k -- 


» A fe P(u) 


) < o 


( 


if Aku = 


4 - 


* vlfeP(u) 


> and 




Aku = 


a k = 


* A fc P(u) 


) > o 


1 


if Aku = 


4 = 


» AfeP(u) 


> and 




A k u = 


a k -- 


=> A k F(u) 


) < o 
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Theorem 2.14 Let m,M G N with to < M . Assume N C W C X m and the 

sequence of pairs {aj G R | < a£, k G N} are topologically self- consistent. 
Fix an integer r > to. Let q = J2*kLr+i Qk<j>ki suc -h that q G n^L J , +1 [a s T, ajj"] and 
qk = for k > M. Let 

N := N x [a m+1 , a+ +1 ] x [a m+2 , a+ +2 ] x ■ • • x [a~, o+] 



Consider the dynamical system induced by §V\ 
If for some j G {to + 1, . . . , r}, dir(j) = 0, then 

CH„(hxv(N)) = 0. 

Assume that for all j G {m + 1, . . . , r}, dir(j) ^ 0, and let d be the number of 
j G {to + 1, . . . , r} such that dir(j) = 1, then 

CH s+d (Inv (N)) = CH s (Inv (N)). (18) 



Proof. By Theorem 2.12, N is an isolating block. 

We will present the proof of the second part of the theorem, only. 

Assume that for all j G {to + 1, . . . , r}, dir(j) ^ 0. Let J := {j | m < j < 
r, dir(j) = 1}. Set 

L-:= (l~x f[ [a^,a+]\u 

\ k=m+l J 

I j-l r 

u [ Nx n [«fe,^]x{a±}x n [ a -, a +] 

j£j \ k=m+l k=j+l 

and 

Z+:= (Vx 

\ k=m+l J 

u ( Arx ii [«*>ofc]*{4} x n 

,30,7 \ fe=m+l fc=J + l 

Let ip:Ex R r — > R r , be any flow generated by 

d k = A k F(u + q) fc = l,...,r 

where a k = A k u and fieWffl IIk=m+i > ffl fc1 ■ 

Clearly, if P r u G Lr , then ^((0, e), P r it) ^ TV for small e > 0. Similarly, if 
P r u G L+, then <p((-e, 0), P r it) £ N for small e > 0. 

Let u — (w,v) G W ffi nl= m +i[ a fc j a fc ]• If w G cW, then either u> is in cW 
or v is in <9n[- =m +i[ a fc j a k ]■ Therefore, dN = L + U L~. 
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Thus, the Conley index of Inv (A) is given by 

CH* (Inv (A) = H*(N,L~). 

A simple argument using the Mayer- Vietoris sequence gives the desired homol- 
ogy groups. I 

Observe that C4a implies that dir(fc) = —1 for all k > m. Therefore, one 
has the following result. 

Corollary 2.15 Assume A C W C X m and the sequence of pairs {aj £ IR | 
a k < a t ' ^ ^ ^1 are topologically self-consistent and satisfy C4a. Fix an 
integer r > m and let 

A := A x [a m+1 , a+ +1 ] x [a m+2 , a+ +2 ] x • • • x [a~, a+]. 

XTien 

CiJ* (Inv (A)) S Off* (Inv (AT)). 

The following theorem is used for all the results described in the Introduc- 
tion. 

Theorem 2.16 Assume A C W C X m and the sequence of pairs {aj: £ M | 
< at , k £ N} are topologically self- consistent and satisfy C4a. Assume 



CHj(lnv (A, tp)) 
for some I, then there exists 



Z if 3=1, 
otherwise, 



u*£Nx H [a kl a+], 

k=m-\-l 

a fixed point for the partial differential equation (Qj. 



Proof. Combining Theorems 2.14 and [2.10| , immediately gives that for each 
r > M there exists a fixed point 



z r £ A x Yl K > 

k— m+1 



for the Galerkin projection onto the first r coordinates. 

Since A x IlfcLm+ifafc > a k ] ^ s com P ac t the collection {z r \ r = m + 1, m + 
2,...} contains a limit point u* . From the continuity of P n o F on W © 
rifeLm±il a fc i a t ] ^ follows that P„ o F(u*) = for each n £ N. By Propo- 
sition 2.4 u* is a fixed point for (|j|). I 
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2.4 Remarks on Related Work 



We are aware of at least two other results that are closely related to the methods 
described in this Section. The first is work of L. Cesari j^] from the early 60's 
which in spirit is very similar to ours. His method can be characterized as follows 
[ p2[ . Let B be a Banach space. Let X be a finite dimensional subspace of B and 
let P : B — > X be a projection. Let N C B be closed with the property that 
PN = N C X is compact and for every x G N, P^ 1 (x) n iV is closed. Consider 
a continuous map / : N — > £?. The goal is to find fixed points for / by studying 
the behavior of the projection of the map onto X . 

It is obvious that u* is a fixed point of /, if and only if Pu* = Pfu*, and 
u* = Pu* + (/ - P)f{u*). 

Cesari's method applies if and only if the following three conditions are 
satisfied: 

(i) For each x £ N, 

P+(I- P)f : P^{x) n N -> P^(x) n N 
is a contraction. 

(ii) Given condition (i), for each x £ N, there exists a unique u{x) E N 
such that Pu(x) + (I — P)f(u(x)) = u(x). The function u : iV — » TV is 
continuous. 

(iii) There are no fixed points of P/it : A" — * iV on the boundary of N. 

Relating this back to the context of this paper, observe that a fixed point 
for the partial differential equation is a fixed point for any nonzero constant 
time map of the corresponding semi-flow, (iii) is closely related to the condition 
C5. As stated (ii) is not well defined unless (i) holds. C4 is the analogous 
assumption to (i) and differs in two significant ways. A necessary condition to 
have a contraction, is for the stronger assumption of C4a to hold. However, C4a 
is not sufficient. An important point is that we do not make any assumptions on 
the direction of the vector field within N. Thus, condition C4a is in principle 
easier to verify than (i). On the other hand, this makes it clear that we cannot 
guarantee uniqueness of the fixed point given our assumptions. 

The other work is due to C. Conley and P. Fife || and is closely related to 
Theorem p,14|. Formulas of the form ( |l8| ) are classical in the context of product 
systems (see j^J). In || one finds a similar formula, but in that context at the 
parameter value for which one computes the index in the lower dimensional 
system, there is no higher dimensional dynamics defined. However, the higher 
dimensional system is defined for an arbitrarily small perturbation. The key 
idea is that in the proper context the lower dimensional system is normally 
hyperbolic. In this paper we circumvent this type of assumption using isolating 
blocks, C5, and imposing C4. 
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3 Estimates for Kuramoto-Sivashinsky equation 



As the Hilbert space H for the Kuramoto-Sivashinsky equation (Q) we choose 
the subspace of L 2 (— it, ir) consisting of 27r-periodic and odd functions. 
Since u(t, x) is odd its Fourier expansion takes the form 

k— oo 

u(t,x) — bk(t) exp(ifcir) (19) 

k— — oo 

Since u(t,x) is real, b k = b-k- Substituting ( |l9| ) into (]!]) gives the following 
equations 

m— oo 

b k = (fc 2 - Vk 4 )b k + ik b ™h-m (20) 



m— — oo 



Since we are interested in solutions with odd symmetry it follows that bk are 
pure imaginary. Let 

a k := ^/~^lb k . 

Then, a k — — a-fc and ao = which results in the following infinite system of 
ordinary differential equations 

k— 1 oo 

at = fc 2 (l - vk 2 )ak - k ^ a n a,k- n + 2fcJ^ a n a n+ k k= 1,2,3,... (21) 

n— 1 n— 1 

We will use these equations to draw rigorous conclusions about Kuramoto- 
Sivashinsky by finding self-consistent apriori boun ds for (Q) that satisfy the 
stronger condition C4a and then applying Theorem 2.16| . To do this, however, 
we will need to understand the errors contributed by ignoring the higher modes 
and the errors introduced by the use of floating point arithmetic. 

Let m, M G N be fixed with m < M. Let W C R m and {a^ G K | fc S N} 
satisfy conditions CI - C3 with the added constraints that 



fe=i 

and 

o*=±§. fc>M (22) 

for some constant C s > and integer s > 1. 

Though technically incorrect, it is perhaps useful for the reader to think of 
the numerical approximation of the dynamics being computed with respect to 
the finite dimensional system 

fc-l M-k 
hk = fc 2 (l - vk 2 )a k - fc ^2 a n0.k-n + 2fc ^2 «n«ii+fc k — 1, . . . , m (23) 

n— 1 n— 1 
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where a& = (a^ + a + )/2 for k = m + 1, . . . M. In doing so it becomes clear 
that there are essentially three levels of approximation that need to be dealt 
with. The first involves the terms in the infinite tail {a k | k > M}. These are 
completely absent from ( |23| ) and therefore must be absorbed as a fixed error 
term (think of the term E{x) in (p"l])). The power decay rule ( |22| ) will be used 
to determine this quantity. The second, involves the terms {a k \ m < k < M}. 
In principle, one could set m = M, however our strategy is to try to obtain 
better estimates for these terms than can be expected by the general decay 
of (B2). However, these terms act as constants and hence can be viewed as 



parameters for the system (23). Finally, the terms {a k \ k = 1, . . . , m} are the 
actual variables for the dynamical system being studied. It should also be kept 
in mind that we need to lift the index information, and therefore need to be 
able to verify C4a for all k. 

Of course, our goal is that of rigorous computations. Therefore each of the 
above mentioned ai is actually an interval. The intervals associated with {a k \ 
k = 1, . . . , m} are essentially determined by the floating point approximations. 
For k > m, the intervals are 

a k = [ak,afc]- 

We will let 

\a k \ := max{| as :|,|a + |}. 

To compute the above mentioned errors we return to ( pl| ) and observe that 
in addition to the linear part there is a finite sum of terms 

fe-i 

FS(k) = £ (24) 



and an infinite sum of terms 



n=l 



oo 



IS{k) = ^ (25) 

71=1 



Obviously bounds on these terms are necessary. 
3.1 1 < k < M 

Since FS(k) is a finite sum and we have already chosen the interval values for a n , 
we can explicitly compute FS(k). Perhaps it is worth noting that to evaluate 
FS(k) only involves the intervals {a n \ n — 1, . . . , M — 1} which are chosen in 
such a way that we expect them to be reasonably good approximations of the 
actual terms. 

Lemma 3.1 Assume 1 < k < M. Then, 

M-k M | 



IS(k) C y^a n a k+n + C s >^ -^-[-1,1] 

fc+i ^ ' 



n=l n=M-k+l 



(k + M+ l) s {s- l)M 



r 2 

-hM] 
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Proof. By definition, 

M-k M oo 

is(k) = 

n=l n=Af-fc+l n=M+l 

With regard to the second sum 

M M ^ 

E «» a fc+» c E K U + s n)s 

Finally, the third sum produces 

E a„a fc+n c £ -£^-^[-1,1] 

n=M+l n=A/+l v ; 



r 2 



(fc + M + l) sL 4^ ™ s 

v ; n=M+l 

c — f-i 11 



In above derivation we used the following estimate 

dx 1 



E - 



< 



n s J M n s (s-UM 3 - 1 



Remark 3.2 This estimate and some of those that follow can be improved by 
noting that 



1 f°° dx 

E — — ^< 



„-w +1 n s {n + k) s J M x s (x + k) s ' 

Of course, the right hand side has an explicit rational expression, but it is rather 
complicated for large s and so was not utilized here. 



A simple extension of Lemma 3.1 leads to the following corollary 



Corollary 3.3 Let 1 < k < m. Then, 

oo M-k M I I 

y~] a n a n+k C Y a n a n+k + C s S~] " [-1, 1] + 

z — ' * — ' f-^ (k + n) s 

n—m—k+1 n—m—k+1 n=M— fe+1 

r 2 

— hi, 11 

(fc + M+l) s (s-l)M s - lL ' 1 



Observe that collorary 3.3 estimates the error in the vector held due to the 



Galerkin projection, namely 

00 

Akip + q)~ A k F(p) = 2k Y a na n +k (26) 



n— m— fc+1 
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3.2 k > M 

Throughout this section it is assumed that k > M. Let 

e(k) •= i 1 if k is even ' 
V ; ' I if fc is odd. 

Lemma 3.4 Let M < k < 2M . Then, 

[k/2] k-M-1 | | 

FS(k)c2 

e(k)a 2 k/2 +2C s TT^lt- 1 ' 1 ]- 

n=fc-M n=l * ' 



Proof. Expanding (|24|) gives 

fe-AT-l Af fe-i 

Qn&k—n "T / &nQ>k—n "T / &n&k—n 
n=l n =k-M n=M+l 

k-M-1 | | M 

c 2C7 s E TFT^t- 1 ' 1 ^ E 

»=1 V ' n=k-M 

[_k/2\ k-M-1 | | 

C ^ a n a fe _„ + e(fc)4 /2 +2C s £ 

n=k-M n=l \ ' 



Lemma 3.5 Let k > 2M . Then, 

C ( 2 S+1 (7 4 s C 2 s \ 

c (^^t g + (2M ; 1)S+1 + (s _ s 1)MS J hi. n 

Proof. From (§|) it follows that 

fc-i 



FS(k) = ^a„a fc _„ 



n=l 
Lfe/2j 

= 2 a n ak- n + e(k)al /2 

n=l 

M |fc/2j 

= 2^a„a fc _„ + 2 2J UnO-k-n + e (fc) a fe/2- 
n=l n=M+l 

Each of these terms will be estimated separately. The first one results in: 

M M 



E a » afe -« c E (fc_MV hl ' 1] 
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c M 



2 S C, W 



On | 
ra=l 

,1/ 



c ^WTT) hl ' 1] ^ KI 



The second term leads to 



[k/2] [k/2] l 

n=M+l n=M+l v ' 

r 2 Lfc/2J - 

= k± y t [-1 11 

n=M+l v ' ; 



civ 1 21 

k s 

C 2 2 s f°° dx 



^ — ^ n s 



k s 1 — ' n 

n=M+l 



C 



C k s J M x^ ^ 

= bvf [-1 11 

k s (s- l)M s ~ l[ ' J 

/-(2os-l 

C s f-1 11 

k s (s-l)M sl ' J 



Finally, the third term gives rise to 

^-[-1 

/fc 2s L ' 1 fc"- 1 (2M+ l) s+1 



^2o2s i C 2 A S 

e(k)al /2 C -^-[-1, 1] C - [-1, l]. 



Turning now to the infinite sum we can obtain the following estimate. 
Lemma 3.6 Let k > M. Then, 

C fe"-i(M + 1) ((M + 
Proof. From (§|) it follows that 

M oo 



M \ 

EKI H,i] 

71=1 / 



is(k) = 



ri=M+l 
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As in the previous case, each term is treated separately. 

n=l n=l ^ ' 

„ AT | 

Gc ^ \a r 



V — [-1 H 

k s ^ (l + n/ky [ ' J 

n— 1 v ' ' 

c M 

V ' 71 — 1 



The remaining term leads to 



V a n a k+n = C 2 S V — 

z — - t m s (k + m) s 



n=M+l n=M+l 



c cL^r * ,_ u , 



(M + 1) S J M (fc + a:)* 

g, 2 

(M + l) s (s- l)(fc + M) 
C 2 



c — [-1 11 

(M + l) s (s-l)fc s - l[ ' J ' 



3.3 Refining the Self-Consistent Bounds 

The proof of our results obviously depends on having good self-consistent bounds 
and the precision of the final result is determined directly by these bounds. For 
this reason it is important to have a process by which these bounds can be 
improved. With this in mind consider an initial sequence of bounds {a^} which 
defines the sets 

TCI OO 

W=Y[[a^,at] and V= JJ [a^,a+]. 

k—1 k— m+1 

We will also assume that 

1 < vrr? . 

This condition means that the Fourier modes for k > m are linearly stable. 

We shall describe the refinement procedure under the assumption of C4a. 
In particular, we need that our sequence {a^ } satisfy dir(fc) = — 1 for all k > m. 
We also assume that since we can numerically solve (|4|), that the estimates for 
W are reasonably good. 

We will inductively adjust Ofc±, for fc = m+1, . . . , M, beginning with a* +1 , 
as follows. Let a <G W V such that a& = at. To satisfy C4a requires that 
d fc < 0, i.e. 

fc 2 (l - vk 2 )a+ - kFS(k) + 2kIS(k) < 0. 
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This is equivalent to requiring 



H > 



2IS(k) - FS(k) 



(27) 



Of course, our goal is to make as small as possible within the constraints 
imposed by the approximations. Since we are iteratively improving our bounds, 
it is reasonable to assume that a worst case equality is the best guess at this 
stage in the procedure. Note, we are not claiming a proof at this point, we 
just are seeking good bounds which later will be verified to be self-consistent 



bounds. So using Lemma 3.1, define fjr to be bounds for 2IS(k) — FS(k), 



M-k M | | 

[/^ + ]:=2£a„a fe+n + 2C s £ -™ 

n= M-k+l V ' 



■1,1] 



71=1 



2C 2 



(k + M + l) s (s- l)M 



fc-i 

Pit - 1j 1] ~ ^ «n«i-n- 



The new value of is given by 



ft 



k ' k*{v-k- 2 y 

A similar argument suggests setting 



fe ' k*{ v -k- 2 y 

This approach works up to k = M. Recall that for k > M, we set a J = 
±C s /k s . Here our goal is to improve the power of convergence, i.e. we want 
to increase s. Again, since we are trying to satisfy C4a the basic inequality 
which needs to be satisfied is (p7f). The estimates for FS(k) and IS(k) for 
k > M obviously are crucial here. However, we had two sets of estimates one 
for M < k < 2M and the other for k > 2M. Thus, we need to choose the worst 
of both estimates. This is done as follows. 

Given an interval J C M let 

1 7" | := sup|x|. 

x£l 



With the estimate from Lemma p.4| in mind define 

Di(k) 



[k/2] 



2 ^2 a nak-n + e(k)a 



2 

k/2 



n —k-M 



k-M-l 



2C S J2 



^ (k- n) s ' 

n=l y ' 



Combining this with the estimate on IS(k) given by Lemma and multiplying 
by k s ~ l leads to the following definition 



2C S 



M 



(M + 1) V (Af+lJ'-H*-!) 



+ y \a n \ ) + max k s ~ 1 D 1 (k). 

l) 5_1 (s-l) ^ M<k<2M ^ ' 
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Turning now to the bounds for k > 2M, Lemmas 3.5, |3.6| , and again multi- 
plying by fc s_1 suggests setting 



Do 



2Ce / 2C S 



(M + i) y(M + iy-i( s -i) _ t 



I2>' 



„s-l 



2M + 1 ^ 1 1 (2M + l) s +! (s - 1)M S 



From Lemmas 3.4, 3.5, and 3.6 we obtain the following result. 



Corollary 3.7 For k > M and D s := max{L»i, D 2 } 

\-FS(k) + 2IS(k)\<^ I . 

We will use this corollary to update the decay rate for the tail terms. Again, 
we want (^) (which gave us C4a) to hold for all k > M. It is sufficient that 

+ \-FS(k) + 2IS(k)\ 
a k > k 3 (v-(M + l)- 2 ) 

and therefore it is sufficient that 

+ D, 11 1 D H 

at > 



k k^ 1 k 3 v-(M+l)- 2 k s + 2 u-(M + l)- 2 

There is a similar inequality for . 

Setting this to an equality we can define 

± , C s +2 , „ D s 



<* l=± W* WherC C '+ a:= ,-(M + l)-» < 28 > 



for k > M. 



4 A Typical Proof 

This section describes the proof of the following result. 
Theorem 4.1 Let 

1 • 1 • „ 
u(x) = — =smx sm2x. 

V ' y/2 8 
For v = 0.75 there exists an equilibrium solution u*(x) to (Qj such that 

| \u* - u\ \ L 2 < 0.052 and | |tt* - u\\ c o < 0.05 
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Stationary Solution to Kuramoto-Shivashinsky v = 0.75 




Figure 2: The computed function u{x) from Theorem 4.1 



The reader should observe that this is a weaker version of Theorem 1.2 



However, we present its proof since it contains all the essential features, but 
with a very low dimensional approximation. The rest of the results described 
in the introduction were proved in a similar manner. 

The first step is to choose m the dimension of the Galerkin approximation 
and M the level to which we make specific choices for the {a^ }. For m = 2, 
23|) reduces to the system of equations 



1 



ai 
a-2 



= -8a, 



- 2a\. 



(29) 



A simple algebraic computation shows that this system has exactly three fixed 
points: (0, 0)-unstable, with one-dimensional unstable manifold and two at- 
tracting fixed points u± — (±-^=,— |). Theorem 4.1 is obtained by studying 

the dynamics of (|29| ) in a neighborhood of (^=, — |). 

The next step is to obtain self-consistent apriori bounds for (|2l]). It is un- 



realistic to expect that goods bounds can be obtained immediately. Thus, we 
make a reasonable guess for bounds and then try to improve them. Let 



W = [ 



1 

71 



0.1, 



1 

71 



0.1] x 



1 1 



0.1]. 
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The initial estimates for [a k , d£] are given in Table |^. The formula used to 
derive these initial estimates will be presented in Section |[ The reason for 
delaying the presentation is to emphasize the fact that the initial estimates are 
only estimates. Obviously, choosing good estimates allows for faster convergence 
and choosing terrible estimates will probably result in a failure of convergence. 



initial 



3 
4 
5 
G 
7 
8 
9 
10 
> 10 



[-0.157542, 0.157542] 
[-0.0463226, 0.0463226] 
[-0.0183725, 0.0183725] 
-0.00871023, 0.00871023] 
-0.00465407, 0.00465407] 
-0.00271036, 0.00271036] 
-0.00168454, 0.00168454] 
-0.00110173, 0.00110173] 
[-1,1] • 10.9915/fc 4 



Table 2: Initial estimates for the intervals [a k , a£]- 

Beginning with the data in Table ^, the refinement procedure described in 
Section |3.3| is used to update for k > 2. After three iterations one obtains 
the estimates given in Table ||. It can now be checked that W and {aj:} from 
Table form self-consistent bounds. 



k 


final [a k , a£] 


3 


[0, 0.021055] 


4 


[-0.00192301, 0] 


5 


[-1.8253 x 10-° 7 , 0.000141734] 


6 


[-9.85549 x 10~ 06 , 8.64999 x lO" 09 ] 


7 


[-6.55526 x 10- 10 , 6.42034 x lO" 07 ] 


8 


[-4.03088 x 10~ 08 , 9.30 9 92 x 10~ n ] 


9 


[-3.51558 x 10- 10 , 2.79203 x lO" 09 ] 


10 


[-1.11597 x 10~ 09 , 9.71368 x 10~ 10 ] 


> 10 


[-1,1] ■ 10285.3/A: 10 



Table 3: Estimates for the intervals [a k , at] representing self-consistent apriori 
bounds 

What should be clear at this point is that the uncertainty contributed by 
the terms not in the Galerkin projection are extremely small. Obviously, at this 
point most of the uncertainty is due to the size of W. 

Having controlled the errors from the Galerkin truncation, the next step is 
to obtain an isolating block N which is topologically self consistent with W 
and {ak}- In determing N we use the vector field (|29|). Of course, the correct 
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equations are given by (pl|): 



a x = -oi + 2y^o n a n+1 



4 

n=l 



2 



Using Corollary 3.3 one obtains the following bounds on the errors, e*, i = 1, 2, 
ei = [-0.00955626, 7.06005 ■ 10~ 10 ] and e 2 = [-1.544 • 10~ 8 , 0.0697171]. 
Thus, the equations for which the isolating neighborhood should be found 

are 

1 , . 

a\ = ~ai + 2aia 2 + t\ (30) 

a 2 = -8a 2 -2al + e 2 - 

The construction of the isolating block around (I/a/2, — 1/8) is easier if one 
works in coordinates determined by the eigenfunctions of the linearized equa- 
tions at the fixed point. The eigenvalues are 

Ai = -4 + 2V3, A 2 = -4-2V3 

with corresponding unit eigenvectors 

vi = , = (1,-=)', «2= = (1,-=)*. 

V15 - 8^3 V2 V15 + 8V3 v2 

Let T be the affine change of variables from (x,yY in this new basis to the 
original variables (ai,et 2 )*. Then, on the set T~ 1 {W) j30| ) becomes 

x = (-A + 2V3)x + f 1 (x,y)+h (31) 
y = (-4-2V3)j/ + / 2 (a;,J/) + e2 

where /i and / 2 are polynomials containing only terms of degree two and e\ , e 2 
are obtained from e±, e 2 by the transformation T. In particular, 

ei = [-0.0110098,0.0152184], e 2 = [-0.0764461,0.00397075] 

Set W = [-0.0748016, 0.0748016] 2 , then T{W) C W. Thus, an isolating 
block N C W (which satisfies the error constraints) will give rise to an isolating 
block N = T{N) C W such that N, W and {at} are topologically consistent. 
Since for each fc, dir(fc) = — 1, 



CH 3 (lnv(N)) 



Z if j = 0, 
otherwise. 
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Thus, by Theorem 2.16, there exists the desired fixed point in 

oo 

u*eNx l[[a^,a+], 

k=2 

Thus, all that remains is to construct N. In the (x,y) coordinates the error 
constraints become 

ft{x, y) + h C (b x ,B a ) := (-0.0334071, 0.0376157) for (x, y) € W 
f 2 (x,y) + e 2 C {by, By) := (-0.0764461,0.00397075) for (x,y) G W 



b x 


B x ~ 


X 






Ai. 


A2 A2 



This implies that bounds on the derivative are given by 

b B b B 

\ 1 (x+^) <x< \ 1 (x+-^) and \ 2 (y + JL) < y < \ 2 (y + -*■). 
Ai Ai A2 A2 

Observe that A,; < 0, hence the box 
N = 

= [ - 0.0623385, 0.0701918] x [-0.0132425, 0.00353264] c W 
is an isolating block. 

5 Obtaining the Initial Estimates 

Before beginning this section we want to once again emphasize that the proofs 
of the theorems in this paper are in principle independent of this section. On the 
other hand, good initial guesses greatly improve the speed of convergence. The 
estimates described in what follows apparently provide excellent initial values 
for the self-consistent bounds. 

We will follow the arguments from |l(| to produce estimates for errors in the 
Galcrkin projection. Kuramoto-Sivashinsky can be written in the form 



vAu - A 1/2 u + 2B(u, u)=0 



where 



A-— A 1 ' 2 --— B(uv)-u— 



While A 1 / 4 ^ it is still the case that 

' ox ' 

l^„b=lgl, 

To simplify the notation, let 

\u\ = \u\ 2 , \\u\\ = lA^uU 
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Since we are interested in bounded invariant sets we can without loss of 
generality assume the following apriori bounds for the invariant set under con- 
sideration: 

|«(*)|<P0, K*)||<Pi, fort>T(«(0)) (32) 
We will make use of the following inequality jl(j, Lemma 1.4] 

\(B(u,v),w)\ < v^|«| 1/a ||«ir /a ||«|||»u| (33) 

The eigenvalues of A are A„ = n 4 , n G N and the corresponding complete 
family of orthonormal eigenfunctions are {-^ sm(nx)}. Let P = P m be an 
orthogonal projection on first m eigenfunctions and set Q = I — P. 

Using the decomposition 

p = Pu, q = Qu 
and the same abuse of notation the equation for q is 

q = -vAq + A 1/2 q- 2QB(u,u). (34) 



Theorem 5.1 Under the assumptions stated above if m is large enough such 
that A m +i > ~^ j then 

limsup|g(t)| < _i . 

W("-*m+i) 

Proof: Beginning with (|m|) and taking a scalar product with q gives 
(^\q) = ~v{Aq\q) + (A^ 2 q\q) - 2(B(u,u)\q). 

Therefore, 

lj t \q\ 2 < -v{Aq\q) + (A 1/2 q\q) + 2\(B(u,u)\q)\ (35) 

Observe that 

(Aq\q) = (A^ 2 q\A^ 2 q) = {A^qf (36) 

and 



(A^ 2 q\q) = J2 

n— m+1 

- x~ 1/2 \ 1/2 xV2u 

— A m+1 2-^i m+l A n \1n 



n— m+1 

n=m-\-l 
oo 

^ A m+1 /\2 A »l 9 "l 
n— m+1 



KS\A l/2 q \ 2 07) 
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Thus, (f35|) (|36|) and (|37|) imply that 

j t \q\ 2 <-2v\A^ q f + 2\;^ 1/2 <l\ 2 +M(B(u,u)\ q )\. (38) 



From (|33D it follows that 

\(B(u,u), q )\ < ^m^imi 3 / 2 !?!. 



While, (|32|) implies that for t > T{u) 

\(B{u,u) iq )\<V2pl /2 p\ /2 \q\ 



Therefore, ( p8| ) becomes 

~M a < -2(*- A^ 2 )!^ 2 ,! 2 + 4V2^ 2 p?/ 2 | 9 | 

Observe that l^l 1 / 2 ^ 2 > A m+1 |<7| 2 and by assumption v — A m ^ 2 > 0. Thus, 
= {^pTpT 2(1/ - \-^ll)\ m ^\ q \)\ q \ 



Thus, for 

1,1 > 4 ^' 2 

^|q| 2 < and hence, 



2 ( i/ - A m ^ 2 )A m+ i 



. /X 1/2 3/2 

limsup| g (t)|< /^V; • 



Remark 5.2 Because an orthonormal collection eigenvectors were used for the 
calculations in this section, the coefficients qk and au differ by a scaling, i.e. 

qu = -2s/7Ta n . 



Theorem 5.1 can be used as follows. For a fixed v, numerical experiments 



can suggest values for po and p\. For m < k < M one can use the formula 



. J , , Pi , 4^2^/90^1 
a fc :=mm ±p 0) ± ¥ ,± fc> _ fc _ 2) 



For k > M, one defines 



04 - 



i/- (M + l)- 
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